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Abstract 

We calculate the partonic structure of a constituent quark in the Next-to-Leading 
Order framework. The structure of any hadron can be obtained thereafter using a 
convolution method. Such a procedure is used to generate the structure function of 
proton and pion in NLO, neglecting certain corrections to Aqcd- It is shown that 
while the constituent quark structure is generated purely perturbatively and accounts 
for the most part of the hadronic structure, there is a few percent contributions coming 
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from the nonperturbative sector in the hadronic structure. This contribution plays 
the key role in explaining the SU (2) symmetry breaking of the nucleon sea and the 
observed violation of Gottfried sum rule. These effects are calculated. We obtained an 
Excellent agreement with the experimental data in a wide range of x = [10 -6 , 1] and 
Q 2 = [0.5, 5000] GeV 2 for the proton structure function. We have also calculated Pion 
structure and compared it with the existing data. Again, the model calculations agree 
rather well with the data from experiment. PACS Numbers: 13.60 Hb, 12.39.-x, 
13.88 +e, 12.20.Fv 

1 INTRODUCTION 

Our knowledge of hadronic structure is based on the hadron spectroscopy and the Deep 
Inelastic Scattering (DIS) data. In the former picture quarks are massive particles and 
their bound states describe the static properties of hadrons. On the other hand, the in- 
terpretation of DIS data relies upon the quarks of QCD Lagrangian with very small mass. 
The hadronic structure in this picture is intimately connected with the presence of a large 
number of partons (quarks and gluons). Mass is not the only difference between these 
two types of quarks; they also differ in other important properties. For example, the color 
charge of quark field in QCD Lagrangian is ill-defined and is not gauge invariant, reflecting 
the color of gluons in an interacting theory; whereas, color associated with a constituent 
quark (CQ) is a well defined entity. It is shown that [1][2] one can perturbatively dress 
a QCD Lagrangian field to all orders and construct a CQ in conformity with the color 
confinement. From this point of view a CQ is defined as a quasi-particle emerging from 
the dressing of valence quark with gluons and q — q pairs in QCD. 

Of course, the concept of CQ as an intermediate step between the quarks of QCD La- 
grangian and hadrons is not new. In fact, Altarelli and Cabibo [3] used this concept in the 
context of SU(Q)XO(3). R.C. Hwa in his work used the term valon,for them, elaborated 
on the concept and showed its applications to many physical processes [4]. In Ref.[2] it 
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is suggested that the concept of dressed quark and gluon might be useful in the area of 
jet physics and heavy quark effective theory. Despite the ever presence of CQ no one has 
calculated its content and partonic structure without resorting to hadronic data and the 
process of deconvolution. Thus, the purpose of this paper is threefold: (a) to evaluate 
the structure of a CQ in the Next-to-Leading Order (NLO) framework of QCD, (b) to 
verify its conformity with the structure function data of nucleon and pion for which ample 
data are available, (c) to include additional refinments that are needed to account for the 
violation of Gottfried Sum Rule (GSR) and the binding effect of constituent quarks to 
form a physical hadron. 

2 FORMALISM 

2.1 Moments of Parton Distribution in a CQ 

In this subsection we will utilize the extended work done on the development of the NLO 
calculations of moments, to calculate the structure of a constituent quark. What will follow 
in the rest of this subsection is not a new next-to-leading order calculation. However, we 
find it interesting to explore the existing calculations in the valon framework. 
By definition, a CQ is a universal building block for every hadron; that is, the structure 
of a CQ is common to all hadrons and is generated perturbatively. Once its structure is 
calculated, it would be possible, in principle, to calculate the structure of any hadron. In 
doing so, we will follow the philosophy that in a DIS experiment at high enough Q 2 , it is 
the structure of a CQ which is being probed, and that at sufficiently low value of Q 2 this 
structure cannot be resolved. Thus, a CQ behaves as a valence quark, and a hadron is 
viewed as the bound state of its constituent quarks. Under these criteria, partons of DIS 
experiments are components of CQ. The structure function of a U-type CQ at high Q 2 
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can be written as [4]: 

Q 2 ) = \z(G± + Gu) + \z{Ga + G A +G f _+G,) + ... (1) 

y u u y u u u u 

where all the functions on the right-hand side are the probability functions for quarks 
having momentum fraction z of a U-type CQ at Q 2 . A similar expression can be written for 
a D-type CQ. Following Ref. [4], we define the singlet (S) and nonsinglet (NS) components 

of the CQ structure functions as: 

/ 

° S = E( G ,% +G«) = G f + (2f- l)G uf (2) 

CQ CQ 

1=1 

f 

g ns = J^iGn - G *) = Gf — G uf (3) 

CQ 

1=1 

where Gf is the favored distribution describing the structure function of a quark within 
a CQ of the same flavor, while the unfavored distribution, G u f, describes the structure 
function of any quark of different flavor within the CQ. We have used, /, for the number 
of active flavors. Solving (2) and (3) for Gf and G u f, we get, 

G f = ±(G s + (2f-l)G NS ), (4) 

G uf = ±(G s -G» s ). (5) 

Having expressed all the structure functions of a CQ in terms of Gf and G u f, we now go 
to the N- moment space and define the moments of these distributions as: 

M 2 (N,Q 2 )= f 1 x N - 2 F 2 (x,Q 2 )dx (6) 
Jo 

Mi(N,Q 2 )= f 1 x N - 1 G i {x,Q 2 )dx (7) 
J o 

where, the subscription i stands for S or NS. Due to charge symmetry, in the following, 
we only write the CQ distribution for proton. 

In the NLO approximation the dependence of the running coupling constant, a s , on Q 2 is 
given by: 

<*s(Q )«— — qT~), (8) 



with [3q = | (33 — 2f) and (3\ = 102 — ^M. The evolution of parton distributions is 
controlled by the anomalous dimensions: 

where i,j = q,g; and ^/( ' 1 ) 7V are one and two loop anomalous dimension matrices. The 
nonsinglet moments in the NLO are as follows: 



47T 2/^o 2/^ asCQo) 

The anomalous dimension matrices, 7(°' 1 ) Ar ,are given by, 



(9) 



(0,l)iV (0,1)JV 

7 (o,Div = I ^ i (10) 

' (0,1)JV (0,1)7V « 

igq 199 



These matrices, govern the moments of the singlet quark and sector as given by: 
M S (N,Q>)\_ o^m ^ N 1 as(Qt)-a s m N N N 
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M G (N, Q 

where 7^ = 7^)^ — 4^7^°^ and 7^°)^, 7^)^ are one and two-loop anomalous dimension 

matrices, respectively. 1 is the unit matrix. Notice that the last term in Eq. (11), Q+< ►_, 

means that all subscripts are exchanged. The leading order behavior is obtained from the 
first term in the square brackets. A± denote the eigenvalues of the one-loop anomalous 
dimension matrix, j(°) N : 

a£ = \ hT + if 9 )N ± ^ifr-i { T? + ^9 )N ifq )N ] (12) 

These anomalous dimensions and 7^'°^ are calculated up to second order in a s and can 
be found in [5] [6]. p± are as follows: 

p N = ±i ^)N_ x N )/{x N_ x N ) (13) 

In Ref. [4], d±^ is used to denote the leading order anomalous dimensions; these are related 
to A± by: 4 = ^. 



In evaluating these moments, we have taken Qq = 0.283 GeV 2 and A = 0.22 GeV, as our 
initial scales. It seems that evolution of parton distributions from such a low value of Qq 
is not justified theoretically; we shall take up this issue shortly in subsection 2.3. 



2.2 Parton Densities In CQ 

In the previous subsection we outlined how to calculate the moments. In this subsection 
we will give the corresponding parton distributions in a CQ. 

The moments of the CQ structure function, F 2 (z, Q 2 ), are expressed completely in terms 
of Q 2 through the evolution parameter t: 

ln QS 

t = In^r (14) 
In ft 

The moments of valence and sea quarks in a CQ are, 

M valence = M NS (N, Q 2 ) (15) 

CQ 

M^ = jj(M s -M NS ) (16) 



where M ' are given by eqs. (9) and (11). It is now straightforward, to evaluate M valence 

^2 



and M^a at any Q or t. These are calculated numerically for a range of N at fixed Q or 
t. In Figure (1) we present these moments both in the leading and next-to-leading orders. 
For every value of t we fit the moments by a beta function, which in turn is the moment 
of the parton distribution as given below, 

zqvaL(z,Q 2 ) =az b (l-z) c (17) 

CQ 

zq^iz, Q 2 ) = azP(l - Z y[l + r,z + ^°" 5 ] (18) 

The parameters a, b, c, a, j3, 7, 77, and £ are functions of Q 2 through the evolution parameter 
t. The same form as in Eq. (18) is obtained for the gluon distribution in a CQ, only with 
different parameter values. The functional forms of these parameters are polynomials of 
order three in t and are given in the appendix. The goodness of these we fits are checked 



by x 2 minimization procedure. We find that, for all parameters a, b, c, and d, of the 
valence quark sector x 2 /DOF = 0.99, with the standard error of order 1CT 3 . For the sea 
quark and gluon sector, x 2 /DOF carries between 0.91 and 0.98. Of course, the forms (17) 
and (18) are not unique; rather, they are the most simple and commonly used forms. We 
note that the sum rule reflecting the fact that each CQ contains only one valence quark 
is satisfied at all values of Q 2 : 



Jo eg 

Substituting these results into Eq. (1) completes the evaluation of CQ structure function 
in NLO. In Figure (2), we plot various parton distributions inside a CQ. 

2.3 discussion 

In subsection 2.1 we indicated that our initial scales are Ql = 0.283 GeV 2 and A = 0.22 
GeV. The above value of Qq corresponds to a distance of 0.36 fm which is roughly equal 
to or slightly less than the radius of a CQ. It may be objected that such distances are 
probably too large for a meaningful pure perturbative treatment. We note that F 2 \z, Q 2 ) 
has the property that it becomes 5(z — 1) as Q 2 is extrapolated to Qq, which is beyond 
the region of validity. This mathematical boundary condition signifies that the internal 
structure of a CQ cannot be resolved at Qq in the NLO approximation. Consequently, 
when this property is applied to Eq.(20), below, the structure function of the nucleon 
becomes directly related to xGcq_(x) at those values of Qq; that is, Qo is the leading 



order effective value at which the hadron can be considered as consisting of only three 
(two) CQ, for baryons (mesons). We have checked that when Q 2 approaches to Qq the 
quark moments approach to unity and the moments of gluon approach to zero. In fact for 
Q 2 = 0.2839 gluon moments are at the order of 10~ 4 and singlet and nonsinglet moments 
are 0.9992. In fact our results are only meaningful for Qq > 0.4 GeV 2 . 
^From the theoretical standpoint, both A and Qq depend on the order of the moments, N; 
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(19) 



but here, we have assumed that they are independent of N. In this way, we have introduced 
some degree of approximation to the Q 2 evolution of the valence and sea quarks. However, 
on one hand there are other contributions like target-mass effects, which add uncertainties 
to the theoretical predictions of perturbative QCD, while on the other hand since we are 
dealing with the CQ, there is no experimental data to invalidate an N independent A 
assumption. 



3.1 Constituent Quark Distribution in hadron 

So far, the structure of a CQ is calculated in the NLO framework. In this section we will 
use the convolution method to calculate the structure functions of proton, F% (x, Q 2 ), and 
pion. Let us denote the structure function of a CQ by F 2 (z, Q 2 ) and the probability of 
finding a CQ carrying momentum fraction y of the hadron by Gcq_ (y). The corresponding 

h 

structure function of a hadron is obtained by convolution of F^iz, Q 2 ) and Gcq (y), 



where summation runs over the number of CQ's in a particular hadron. We note that 
Gco_(y) is independent of the nature of the probe, and its Q 2 value. In effect Gcq (y) 

h h 

describes the wave function of hadron in CQ representation, containing all the complica- 
tions due to confinement. ^From the theoretical point of view, this function cannot be 
evaluated accurately. To facilitate phenomenological analysis, following [4], we assume a 
simple form for the exclusive CQ distribution in proton and pion as follows, 



3 HADRONIC STRUCTURE 



h 




(20) 



GuuD/p(yi, 2/2,2/3) = %i2/2) m 2/3 <%i + 2/2 + 2/3-1) 



(21) 



GuD/n- (yi, 2/2) = qyivZSivi + 2/2-1) 



(22) 
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where I and q are normalization parameters. After integrating over unwanted momenta, 
we can arrive at the inclusive distribution of individual CQ: 

G ^>= B(„+l!2„, + 2/ ' (1 -^ < 24 > 



1 

B(n+T,u+ 1) 



G o/«- (v) = p ,„ i , „ i ^ (i - v) w - (25) 



Similar expression for G D /^- is obtained with the interchange of ji v. In the above 
equations, B(i,j) is the Euler Beta function. The arguments of this function, as well as, 
I and q of Eqs.(21) and (22), are fixed by the number and momentum sum rules: 

f 1 GcQ(y)dy = l, V f 1 yGcQ(y)dy = l, (26) 

JO h CQ Jo h 

where CQ = U,D,U and h = p,ir~. Numerical values are: /j, = 0.01, v = 0.06, m = 0.65 
and n = 0.35. These parameters are independent of Q 2 and are given in [4] and [7] for 
proton and pion, respectively. In [8], a new set of values for m and n are suggested which 
differ significantly from the values quoted above. Reference [8] uses the CTEQ parton 
distributions in the Next-to-Leading order to fix these parameters. The authors of Ref. 
[8] have found yet another set of values in [9] , which differs from those given in [4] and [8] . 
They attribute these differences to the different theoretical assumptions. We have tried 
a range of values for m, n, v and fi, and checked that the sum rules (26) are satisfied. 
Because of the Beta function in the denominator of Eqs. (23-25) the number sum rule does 
not change by changing m and n and we found little sensitivity on the momentum sum 
rule. In Fact, we have varied both m and n, randomly in the interval [0.2,1.95], and 
checked it against data. We found very little sensitivity to the choices made. For 
the data that we shall analyze, these two parameter formulas will prove to be adequate. 
The differences between parameters of Ref. [8] and Ref. [4] (which we have adopted) can 
be attributed to the different methods of extracting the parameters. We have calculated 
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the parton distributions from the perturbative QCD, whereas Ref.[8] have used CTEQ's 
global fit to extract parton content of nucleon. 

The considerations on the form of CQ distribution that are used here, does not exclude 
other possibilities. However, in our model the quark distribution, say, in a proton is de- 
scribed as the convolution of the CQ distribution, Gcq_, and the quark distribution in the 
CQ. In the moment form the latter two become the product of their moments. Mathemat- 
ically, these two moments can be modified without changing the product. That would, of 
course, modify the original Gcq_ and parton distribution in a CQ. But the parton distri- 
bution has a definite meaning. It is the evolution of a quark as Q increases. At Q = Qo> 
parton distribution is a delta function, implying that the CQ has no internal structure 
that can be probed. Our model is to say that the three valence quark at high Q becomes 
three CQ's at low Q without anything else. Such a physical statement of the model does 

not permit to redefine Gcq_ or the parton distribution in a CQ, although mathematical 

p 

criticism still remains. 

Our choice of Gcq_ is sufficient to fit all the data of deep inelastic scattering in the range 

p 

of interest. It does not mean that a more complicated Gcq_ is ruled out. It is a subject 

p 

that can be revisited in the future when problem arise. At this point our formula for the 
CQ distribution is satisfactory. A comment regarding the moments is in order: The defi- 
nition of the moments are not arbitrary if we are to attach physical interpretation to the 
constituent quark distribution. There are only two U and one D in a proton, which impose 
constrains on the first moment. Also, the sum of the CQ momenta should be one, which 
imposes a constraint on the second moment. If one writes the valence quark moments 
as a product of G n and M NS then M NS must have connection with QCD evolution to 
ensure that there is only one valence quark in every CQ and hence, the constraint becomes 
M^ s = 1. With these constrains, one cannot freely redefine moments of CQ distribution, 
M NS and M s . yet, there are infinite number of CQ moments and we have cited only two 
constraints. Does this leave an infinite horde number of undefined quantities? The higher 
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moments of the CQ distribution are affected mainly by the large x behavior of proton 
wave function, or will affect only the even larger x behavior of the parton distribution. 
Until one focuses on that region and finds inconsistency, one chooses the simplest param- 
eterization. There are similar situations both in classical and contemporary physics. The 
question can be turned around and asked why did Newton stop at \ law for gravitational 
force? Could he have excluded the possible existence of ^ terms with n = 3,4, 5, ...? By 
only citing the existence of two constraints, Newton avoided confronting the infinite horde 
of undefined quantities. Our model is inductive physics. Valons or CQ's are not derived 
by deductive logic from the first principles. In modeling, one finds the simplest way to 
capture the essence of the physics involved. The situation is also present in contemporary 
physics: we could have asked why the Cornell potential for the bound-state problem of 
quark and antiquark does not contain all possible power law terms beside the - and the r 
terms. The answer is simple: it is what fits the particle mass spectrum that counts, again 
in the simplest way possible that captures the essence of the problem. 

3.2 Parton distribution in hadrons 

Having specified the parton distributions in the CQ and the constituent quark distribution 
in proton and pion, it is possible to determine various parton distributions in any hadron. 
For proton, we get: 

/l dy x 
-jG u/p (y)q vaL/u (-,Q 2 ) 

(' dy x 
+ / —G D / p (y)q vaL / D (-, Q 2 ) = u vaL/p (x, Q 2 ) + d vaL/p (x, Q 2 ) (27) 
Jx y y 

f dy x f dy x 

q S ea/ P (x,Q 2 ) = 2 / —G u/p (y)q sea/u (-,Q 2 )+ —G D/p (y)q sea/D (-,Q 2 ) (28) 
Jx y y Jx y y 

The above equations represent the contribution of constituent quarks to the sea and the 
valence quark distributions in proton. In Figure (3), the CQ contribution to parton 
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distributions in p and ir~ is shown. Comparing with proton structure function data, shows 
that the results fall short of representing the experimental data by a few percent (3%-8%, 
depending on Q 2 and x). This is not surprising, for it is a well established fact that there 
are soft gluons in the nucleon. We attribute their presence to the fact that constituent 
quarks are not free in a hadron, but interact with each other to form the bound states. 
At each Q 2 value we have computed the soft gluon contribution using direct comparison 
with experimental data. The final result is parameterized as: 

xg soft = R(Q 2 )x-°™(l-x) 2 ™; (29) 

where, R(Q 2 ) = 2.14 - 6.13* + 6.177t 2 - 1.812t 3 . In figure (4), the shape of xg soft is 
presented for a range of Q 2 values. Splitting of these gluons could in turn produce q — q 
pairs which can combine with the original constituent quarks and fluctuate the baryonic 
state to a meson-baryon state. This also will break the SU (2) symmetry of the sea as we 
shall see next. 



3.3 The Role of Soft Gluon 

In our model there is no room in the CQ structure for the breaking of SU(2) symmetry 
of the nucleon sea. But, after creation of q — q pair by the soft gluons, these quarks can 
recombine with CQ to fluctuate into meson- nucleon state which breaks the symmetry of 
the nucleon sea. In what follows we will compute this component, and its contribution to 
the nucleon structure function, and its relevance to the violation of the Gottfried sum rule. 
We will follow the prescription of [10] and [11]. The above proposed model, is similar to 
the effective chiral quark theory, in which the relevant degrees of freedom are constituent 
quarks, gluons and Goldstone bosons. The coupling of pions produces pion cloud of the 
constituent quark, in which we are interested. The point we would like to make is that, 
there are two types of sea quark and gluon that contribute to the nucleon sea measured 
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in DIS. One type is generated from QCD hard bremsstrahlung and gluon splitting. This 
component is discussed in subsection 2.2 and is associated with the internal composition 
of the constituent quark rather than the proton itself. The other contribution is generated 
nonperturbatively via meson-baryon fluctuation. As we will show, according to our results, 
this component provides a consistent framework to understand the flavor asymmetry of 
the nucleon sea. It also compensates for the shortcoming of perturbatively generated 
parton densities in reproducing the experimental data on F^. In order to distinguish 
these partons from those confined inside the CQ, we will term them as inherent partons. 
This component is intimately related to the bound state problem, and hence it has a non- 
perturbative origin. However, for the process of CQ — > CQ + gluon —> q — q, at an initial 
value of Q 2 = 0.65GeV" 2 where a s is still small enough, we will calculate it perturbatively. 
The corresponding splitting functions are as follows, 

P m(z) = li±^, (30) 

P qa (z) = ±(z 2 + (1 + z) 2 ). (31) 
For the joint probability distribution of the process at hand, we get, 

Qinhix, Q 2 ) = q inh (x, Q 2 ) = AT— A- / —Pq 9 {-)g S oft{y)dy (32) 

where qinh. is the inherent parton density. The splitting functions, and the qi n h(x,Q 2 ) = 
(<linh( x , Q 2 )), are those of the leading order, rather than NLO. We do not expect that this 
approximation will affect the whole structure function, as can be seen from Figure (5). 
In the above equation M, is a factor depending on Q 2 , and Gqq is the constituent quark 
distribution in the proton as given previously. The same process, can also be a source of 
SU(2) symmetry breaking of nucleon sea; namely, u sea d sea , and hence, the violation of 
the Gottfried sum rule (GSR). We will take up this issue in the next subsection. 
The process is depicted in Figure (6). Probability to form a meson-baryon state can be 
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written as in Ref. [11], 

Pmb(x)= j 1 ^- f 1 -F(y,z)R(y,z;x) (33) 
Jo y Jo z 

where F(y, z) is the joint probability of finding a CQ with momentum fraction y, and 
an inherent quark or anti-quark of momentum fraction z in the proton. R(y, z; x) is the 
probability of recombining a CQ of momentum y, with an inherent quark of momentum 
z, to form a meson of momentum fraction x in the proton. For a more general case, the 
evaluation of both of these probability functions are discussed in [12]. An earlier and 
pioneering version was proposed in [13]. In the present case, these functions are much 
simpler. Guided by the works of Ref. [11, 12, 13], we can write, 

F(y,z) = nyGcQ(y)zq mh Xz)(l-y-z) s (34) 

V 

R(y,z;x) = py a z b 5(y + z-l) (35) 

Here, we take a = b = 1 reflecting that the two CQ's in meson share its momentum almost 
equally. The exponent 5 is fixed for the n and A ++ states, using the data from E866 
experiment [14], and the mass ratio of A to n. They turn out to be approximately 18 
and 13, respectively. O and p are the normalization constants, also fixed by data. We 
recognize that the discussion leading to Eqs. (34, 35) and fixing the parameters are based 
on a phenomenological ansatz. The ansatz and the parameters of Eqs. (34, 35) seems 
working remarkably well, and further effects that might be of hadronic origin need not be 
larger than a few percent. 

Now it is possible to evaluate % and du quarks, associated with the formation of meson 
statesfll]: 

d M (x,Q 2 ) = C ^[P w+n (y) + ^ +A o(y)K(-,Q 2 ) (36) 
Jx y o y 

1 dv x 

um(x,q 2 ) = - —Pn-A++(.y)M-,Q 2 ) (37) 
£ Jx y y 

where u w and d w are the valence quark densities in the pion at scale Qq. The coefficients ^ 
and | are due to isospin consideration. Using Eqs. (17), (18), and (25), we can calculate 
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various parton distributions in a pion. Those pertinent to Eqs. (35) and (36) are: 

ulai.^Q 2 ) = [ l G D/ ^{y)u vaL/0 { X -,Q 2 ) d -^ 1 (38) 

Jx y y 

Q 2 ) = C Gd/tt- (y)d va i./D(-, Q 2 )—, (39) 
Jx y y 

U V al./U = U V al./U- (40) 

There are some data on the valence structure function of 7r~ [15], which can serve as a 
check for the proposed model. Defining valence structure function of tt~ as 

Kal. = X<al. = xdZal. (41) 

We present the results of our calculation for F£ al in Figure (7), along with the experimental 
data at Q 2 around 6 Gev 2 . 

3.4 Asymmetry of the Nucleon Sea and the Gottfried Sum Rule 

There are experimental evidences [14], [16] that the Gottfried integral, 

S G = / [it (x) - F?(x)]— = - - - / dx[d(x) - u(x)}, (42) 
Jo x 6 3 jo 

is less than | which is the value expected in the simple quark model. The NMC collabora- 
tion [16], result at Q 2 = 4 GeV 2 is, Sq = 0.235 ±0.026, which is significantly smaller than 
|. There are several explanations for this observed violation of the Gottfreid sum rule 
(GSR). Among them are: flavor asymmetry of the nucleon sea, that is, u sea / J sea ,[17] 
[18], isospin symmetry breaking between proton and neutron, and Pauli blocking, among 
others. One of these explanations fits well within our model. 

It was proposed by Eichten, Inchliffe, and Quigg [19], that valence quark fluctuates into a 
quark and a pion. This explanation fits very well in our model. In other words, a nucleon 
can fluctuate into a meson-nucleon state. The idea is appealing to us and in our model 
it can be calculated rather easily. According to our model, after a pair of inherent q — q, 
is created, a u can couple to a D-type CQ to form an intermediate ir~ = Du, while the 
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u quark combines with the other two U-type CQ's to form a A ++ . This is the lowest 
uu fluctuation. Similarly, a dd can fluctuate into the 7r + n state. Since the A ++ state is 
more massive than the n state, the probability of dd fluctuation will dominate over uu 
fluctuation, which naturally leads to an excess of dd pairs over uu in the proton sea. This 
process is depicted in Figure (6). 

To summarize, there are three sources that contribute to the sea partons in the proton: 
Constituent quarks of the proton, the mesonic cloud of the CQ, and the splitting of soft 
gluon. The combined contribution to the ratio of = is as follows: 

,d. _ d M + d inh , + cQ . s 

V _ jproton — _ • \^"J 

U Um + Ui n h.+CQ 

NUSea collaboration at FermiLab, E866 experiment, has published their results [14], for 
the integral of d — u and | at Q = 7.35 GeV. With the procedure described, we have 
calculated these values at the same Q and for the same range of x, that is, x = [0.02, 0.35]. 
The result of the model is, 

r0.345 

dx(d -u) = 0.085, (44) 



'0.02 

to be compared with the experimental value of 0.068 ± 0.0106. For the entire range in 
x, we get: Jq dx(d — u) = 0.103. The corresponding experimentally extrapolated value is 
0.1 ± 0.018, which is in excellent agreement with our calculations. At Q = 7.35 GeV, this 
gives Sq = 0.264. In Figure (8), d{x) — u(x) and | in the proton are shown as a function 
of x at Q = 7.35 GeV along with the measured results. 



3.5 Proton Structure Function F$ 

We are now in a position to present the results for the proton structure function, F% . In 
Eqs. (17,18) we presented the form of parton distributions in a CQ. Using those relations, 
along with the numerical values given in the Appendix, the structure function of the CQ 
is obtained via Eq. (1). The shape of CQ distributions in the proton is provided by Eqs. 
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(23) and (24). Now, all the ingredients are in place to calculate i*f from Eq. (20). In 
Figure (9) the results are shown at many values of Q 2 . As it is evident from Figure (4), 
sole use of CQ structure to reconstruct the F$ data is not enough. The CQ contribution 
fall a few percent short of representing the actual data. However, as mentioned earlier, 
there is an additional contribution from the inherent partons to F^, which is calculated 
in Eq. (31). Adding this component represents the data rather well as can be seen from 
Figure (9). The data points are from [20]. For the purpose of comparing our results with 
other calculations, we have also included in Figure (9), the GRV's NLO results [21], as well 
as the prediction of CTEQ4M [22]. Notice that we have taken the number of active flavors 
to be three for Q 2 < 5GeV 2 and four, otherwise. In Figure (10), the gluon distribution 
predicted by the model is presented along with those from Ref.[16, 20]. 

4 Summary and Conclusion 

In this paper we have used the notion of constituent quark as a well defined entity being 
common to all hadrons. Its structure can be calculated perturbatively to all orders in QCD. 
A CQ receives its own structure by dressing a valence quark with gluon and q — q pairs in 
QCD. We have calculated its structure function in the Next-to-Leading order framework. 
Considering a hadron as the bound states of these constituent quarks, we have used the 
convolution theorem to extract the hadronic structure functions for proton and pion. 
Besides the CQ structure contribution to the hadrons, there is also a nonperturbative 
contribution. This component is small, being only a few percent. But it is crucial in 
providing a framework to understand and explain the violation of the Gottfried sum rule 
and the excess of d over u in the nucleon sea. A mechanism is devised for this purpose and 
necessary calculations are performed. We have presented the results and compared them 
with all the available and relevant data, as well as with the work of others. We found that 
our results are in good agreement with the experimental data. 
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5 APPENDIX 

In this appendix we will give the functional form of parameters of Eqs. (17, 18) in terms of 
the evolution parameter, t. This will completely determines partonic structure of CQ and 
their evolution. The results are valid for three and four flavors, although the flavor number 
is not explicitly present but they have entered in through the calculation of moments. As 
we explained in the text, we have taken the number of flavors to be three for Q 2 < 5GeV 2 
and four for higher Q 2 values. 

I) Valence quark in CQ (Eq. 17): 

a = -0.1512 + 1.785* - 1.145* 2 + 0.2168* 3 
b = 1.460 - 1.137* + 0.471* 2 - 0.089* 3 
c = -1.031 + 1.037* - 0.023* 2 + 0.0075* 3 

II) Sea quark in CQ (Eq. 18): 

a = 0.070 - 0.213* + 0.247* 2 - 0.080* 3 
P = 0.336 - 1.703* + 1.495* 2 - 0.455* 3 
7 = -20.526 + 57.495* - 46.892* 2 + 12.057* 3 
v = 3.187 - 9.141* + 10.000* 2 - 3.306* 3 
£ = -7.914 + 19.177* - 18.023* 2 + 5.279* 3 
M = 1.023 + 0.124* - 2.306* 2 + 1.965* 3 
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Ill) Gluon in CQ (Eq. 18) 

a = 0.826 - 1.643* + 1.856* 2 - 0.564* 3 
= 0.328 - 1.363* + 0.950* 2 - 0.242* 3 
7 = -0.482 + 1.528* - 0.223* 2 - 0.023* 3 
rj = 0.480 - 3.386* + 4.616* 2 - 1.441* 3 
£ = -2.375 + 6.873* - 7.458* 2 + 2.161* 3 
M = 2.247 - 6.903* + 6.879* 2 - 1.876* 3 
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6 Figure Caption 

Figure-l. Moments of partons in a CQ at Q 2 = 8.5,25,120 GeV 2 . The dashed curve 
represents the variation of leading order moments and the solid curve is that of next to 
leading order. 

Figure-2. Parton distributions in a CQ at a typical value of Q 2 = 20 GeV 2 

Figure-3. Parton distributions in proton and ir~ at Q 2 = 20 GeV 2 as a function of x. 

Here only perturbatively generated partons are presented. 

Figure-4 Soft gluon distribution, xg so jt, as a function of x at various Q 2 values. 
Figure-5. The effect of inherent component to a sea parton distribution in proton is 
presented as a function of x. The dashed-dotted line is the perturbatively generated in- 
constituent quark contribution. The solid line represents the sum of the two components. 
Figure-6. Processes responsible to SU{2) symmetry breaking in the nucleon sea and 
violation of Gottfried sum rule. 

Figure-7. Pion valence structure function as a function of x at Q 2 = 5.5 (Gev/c) 2 . Solid 
curve is the result of the model calculations and the data points are from Ref.[17]. 
Figure-8. The ratio | and the difference d — u as a function of x. The solid line in the 
model calculation and the dotted line is the prediction of CTEQ4M [22]. Data are from 
Ref. [18]. 

Figure-9. Proton structure function F% as a function of x calculated using the model and 
compared with the data from Ref. [20] for different Q 2 values. The dot-dot-dashed line is 
the prediction of GRV Ref. [21] and the dashed line is that of CTEQ4M Ref. [22]. 
Figure-10. The gluon distribution in proton as a function of x at Q 2 = 20 GeV 2 . We have 
also shown the prediction of GRV (dashed-dotted line) and CTEQ4M (dashed line). The 
data points are from HI collaboration. 
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